Application of a new fast multipole BEM for simulation of 2D elastic solid with large number of inclusions

A fast multipole method (FMM) is applied for BEM to reduce both the operation and memory requirement in dealing with very large scale problems. In this paper, a new version of fast multipole BEM for 2D elastostatics is presented and used for simulation of 2D elastic solid with a large number of randomly distributed inclusions combined with a similar subregion approach. Generalized minimum residual method (GMRES) is used as an iterative solver to solve the equation system formed by BEM iteratively. The numerical results show that the scheme presented is applicable to certain large scale problems. The project supported by the National Nature Science Foundation of China (10172053) and the Ministry of Education     

适合于求解边界元方程组的GMRES算法的实用化和并行化研究

为了将GMRES算法应用于大型边界元方程组的求解，采用预条件技术和重正交技术相结合的方法实现了该算法的实用化，然后在实用化的基础上针对迭代算法具有良好并行性的特点，研究了该算法在网络机群环境下的并行化技术。数值试验和分析表明所用的这些技术是行之有效的，对于提高求解速度和增大求解问题的规模是有意义的。     

On the block GMRES method with deflated restarting

In this paper, we consider the block-GMRES method with deflated restarting for solving nonsymmetric linear systems with multiple right-hand sides. We modify slightly the restarted block-GMRES method with deflation of eigenvalues proposed by Morgan to obtain a new one. It is shown that the modified method is mathematically equivalent to the Morgan’s original method. However, the error analysis shows the modified version minimizes the numerical errors during a restart and therefore is better suited if the linear systems have to be solved with high precision. Numerical experiments report the effectiveness of the modified method.     

Fully coupled flow-induced vibration of structures under small deformation with GMRES method

Lagrangian-Eulerian formulations based on a generalized variational principle of fluid-solid coupling dynamics are established to describe flow-induced vibration of a structure under small deformation in an incompressible viscous fluid flow. The spatial discretization of the formulations is based on the multi-linear interpolating functions by using the finite element method for both the fluid and solid structures. The generalized trapezoidal rule is used to obtain apparently non-symmetric linear equations in an incremental form for the variables of the flow and vibration. The nonlinear convective term and time factors are contained in the non-symmetric coefficient matrix of the equations. The generalized minimum residual （GMRES） method is used to solve the incremental equations. A new stable algorithm of GMRES-Hughes-Newmark is developed to deal with the flow-induced vibration with dynamical fluid-structure interaction in complex geometries. Good agreement between the simulations and laboratory measurements of the pressure and blade vibration accelerations in a hydro turbine passage was obtained, indicating that the GiViRES-Hughes-Newmark algorithm presented in this paper is suitable for dealing with the flow-induced vibration of structures under small deformation.     

A Fast,Matrix-free Implicit Method for Computing Low Mach Number Flows on Unstructured Grids

A fast, matrix-free implicit method has been developed to solve low Mach number flow problems on unstructured grids. The preconditioned compressible Euler and Navier-Stokes equations are integrated in time using a linearized implicit scheme. A newly developed fast, matrix-free implicit method, GMRES + LU?SGS, is then applied to solve the resultant system of linear equations. A variety of computations has been made for a wide range of flow conditions, for both in viscid and viscous flows, in both 2D and 3D to validate the developed method and to evaluate the effectiveness of the GMRES + LU?SGS method. The numerical results obtained indicate that the use of the GMRES + LU?SGS method leads to a significant increase in performance over the LU?SGS method, while maintaining memory requirements similar to its explicit counterpart. An overall speedup factor from one to more than two order of magnitude for all test cases in comparison with the explicit method is demonstrated.     

二维弹性新型快速多极虚边界元的最小二乘法

在虚边界元最小二乘法的方程求解中采用新型的快速多极展开和广义极小残值法,提出了一种二维弹性新型快速多极虚边界元最小二乘法的求解思想。基于二维弹性问题原有的快速多极虚边界元最小二乘法的展开格式,通过引入对角化的概念,以更新展开传递格式;相对于原有快速多极算法,该方法可进一步提高计算效率且仍能保证具有较高的计算精度。数值算例说明了该方法的可行性、计算效率、计算精度均较高。     

A parallel computational method for simulating two‐phase gel dynamics on a staggered grid

We develop a parallel computational algorithm for simulating models of gel dynamics where the gel is described by two phases, a networked polymer and a fluid solvent. The models consist of transport equations for the two phases, two coupled momentum equations, and a volume‐averaged incompressibility constraint. Multigrid with Vanka‐type box‐relaxation scheme is used as preconditioner for the Krylov subspace solver (GMRES) to solve the momentum and incompressibility equations. Through numerical experiments of a model problem, the efficiency, robustness and scalability of the algorithm are illustrated. Copyright © 2008 John Wiley & Sons, Ltd.     

A MODIFIED GMRES METHOD FOR SOLVING LARGE NONSYMMETRIC LINEAR SYSTEMS

A modified GMRES method is proposed in this paper, the method replaces the approximation xm obtained by the GMRES method with a new approximation xm which is a linear combination of xm and the wasted basis vector vm 1. The residual norm of the new approximation satisfies a small one-dimensional minimization problem. Relationships between the residual norms of xm and xm are given. We show that the resulting m-step modified GMRES method is better than the original m-step GMRES method in theory and is consi...     

Comparison of GMRES and ORTHOMIN for the black oil model on unstructured grids

This paper addresses an application of ORTHOMIN and GMRES to petroleum reservoir simulation using the black oil model on unstructured grids. Comparisons between these two algorithms are presented in terms of storage and total flops per restart step. Numerical results indicate that GMRES is faster than ORTHOMIN for all tested petroleum reservoir problems, particularly for large scale problems. The control volume function approximation method is utilized in the discretization of the governing equations of the black oil model. This method can accurately approximate both the pressure and velocity in the simulation of multiphase flow in porous media, effectively reduce grid orientation effects, and be easily applied to arbitrarily shaped control volumes. It is particularly suitable for hybrid grid reservoir simulation. Copyright © 2006 John Wiley & Sons, Ltd.